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Abstract 

Superbeam experiments can, in principle, achieve impressive sensitivities for CP violation in 
neutrino oscillations for large 6*13. We study how those sensitivities depend on assumptions about 
systematical uncertainties. We focus on the second phase of T2K, the so-called T2HK experiment, 
and we explicitly include a near detector in the analysis. Our main result is that even an idealised 
near detector cannot remove the dependence on systematical uncertainties completely. Thus ad- 
ditional information is required. We identify certain combinations of uncertainties, which are the 
key to improve the sensitivity to CP violation, for example the ratio of electron to muon neutrino 
cross sections and efficiencies. For uncertainties on this ratio larger than 2%, T2HK is systematics 
dominated. We briefly discuss how our results apply to a possible two far detector configuration, 
called T2KK. We do not find a significant advantage with respect to the reduction of systematical 
errors for the measurement of CP violation for this setup. 
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I. INTRODUCTION 



Neutrino oscillation offers a natural, minimalistic framework to account for the observed 
data in a wide variety of experiments. It has been established as the leading mechanism 
for flavour transitions in solar [l|, 0] and atmospheric Q neutrinos. The typical L/E pat- 
tern expected for oscillations begins to emerge from various data sets like the most recent 
KamLAND data 0, \^ or the first generation of long-baseline disappearance experiments 
K2K @, 3 and MINOS [8|. The only experiment which cannot be accounted for in a three 
flavour oscillation framework is the LSND V e appearance signal (||. The LSND result is 
quite difficult to reconcile with the null results of CDHS 10j and Bugey 11], even when 



one allows for one or more sterile neutrinos. Recently, MiniBooNE failed to confirm the 



LSND signal 12J. For an analysis of sterile neutrino solutions to the LSND result in view of 



the MiniBooNE data see [13[. The status of LSND thus remains unclear and therefore we 
will assume that LSND has a non-oscillation explanation. We will consider only oscillations 
between three active flavours. 

Neutrino oscillations require massive neutrinos. This on its own is one of the strongest 
and first indications for physics beyond the Standard Model. Most models for neutrino 
mass generation point to a very high energy scale far beyond the reach of current and future 
accelerator experiments. The neutrino is thus a unique messenger for otherwise in-accessible 
physics. To fully exploit this potential, new, high precision oscillation experiments are nec- 
essary. These experiments will address the size of #13, the neutrino mass hierarchy, leptonic 
CP violation and whether #23 = vr/4. None of the currently running or approved experiments 
has sufficient sensitivity to achieve an accurate measurement of the neutrino mass hierarchy 
or CP violation [14] . The reason is, that in both cases the effects are quite small and subtle. 
There is a plethora of possible technologies to address these questions. They range from 
third generation superbeam experiments to beta beams and neutrino factories. For a recent 
review see 15 



In this paper we focus on third generation superbeam experiments. These experiments 
are based on conventional neutrino beams from 7r-decay. The experiments are 'super' in 
the sense that they will use proton beams of unprecedented strength around 1 — 4 MW and 
detectors with fiducial masses of several 100 kt. This will allow to collect many thousands 
of v e and V e appearance events (assuming sin 2 29 13 = 0.1). Thus statistical errors will be 
at most of percent size and systematical errors will be no longer be negligible. Nearly 
all previously published sensitivity studies use some ad-hoc 1 value of systematics which is 
assumed to be achieved by means of a near detector. This near detector is not specified nor 
included in the calculation. 

In this paper we carefully investigate a large number of possible contributions to the 
systematical error budget in superbeam experiments. To be specific we use the T2K [16 



1 ad-hoc refers to the fact that this value is usually smaller than any value achieved by any previous 
experiment. 
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experiment as example, focusing mainly on the discovery of CP violation in T2K phase II 



17, 



18 



(T2HK). We briefly comment also on T2K phase I, as well as the T2KK proposal 
where half of the T2HK detector is moved to Korea. Our analysis is based on a realistic 
Monte Carlo study of the detector response and we explicitly include a (though somewhat 
idealised) near detector in the simulation. The goal of this paper is to investigate to which 
extent a near detector can contribute to reduce the impact of various systematical error 
sources. We will show that even an idealised near detector on its own cannot reduce the 
impact of all error terms. Thus additional information on fluxes and/or cross sections will 
be required. 

The outline of the paper is as follows. In Sec. [Til we discuss qualitatively the impact 
of systematics on an appearance experiment. Though highly simplified this discussion will 
allow us to understand many features of the numerical calculations. In Sec. II II I we review 
in some detail uncertainties on neutrino cross sections, since they will be crucial in the 
subsequent analysis. In Sec. IIVI we give a brief description of our experiment simulation 
and the various types of systematics we include in our analysis. More technical details 
can be found in appendix [A] Sec. [V] contains the results of our numerical calculations: in 
Sec. IV Al we consider the CP violation sensitivity of T2HK discussing mainly the impact 
of cross section uncertainties, whereas in Sec. IV Bl we point out under which circumstances 
information on fluxes and an improved near detector are useful. In Sec. IV CI we consider 
the determination of #13 in the context of T2K (phase I) and T2HK (phase II), whereas in 
Sec. IV Dl we discuss the systematics impact for T2KK in comparison with T2HK focusing 
again on CP violation. A summary and some speculative remarks on other high precision 
oscillation facilities follow in Sec. I VII In appendix IA II we give technical details on the 
experiment simulation, whereas in appendix I A 21 we describe the statistical analysis including 
the implementation of systematics. In appendix [B] we show how the full ND/FD simulation 
can be approximated by a much simpler FD-only analysis adopting only very few "effective" 
systematics parameters. 



II. QUALITATIVE DISCUSSION 

Before we plunge into a detailed numerical study, we will present here a simple argument, 
why even an idealised near detector is not sufficient for an appearance experiment. For the 
sake of discussion we introduce a couple of simplifications, for example we consider only total 
rates and neglect some background sources, while we do include them in the subsequent 
numerical calculations, see Sec. [IV] and App. [A] for details. 

The total number of v e and events in near detector (ND) and far detector (FD) can 
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be written as 



nd = l^ro^ (1) 
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Here N is the total normalisation (number of target nuclei), a Va is the charged current cross 
section for u a , e Ua is the detection efficiency for u a (assumed to be identical for ND and 
FD), P{yp — > u a ) is the probability for a neutrino of flavour j3 to oscillate into flavour a, 
$> Ug is the initial neutrino flux, and L is the distance from the detector to the source. For 
the v e signal we include here the intrinsic u e beam contamination and the background from 
neutral current (NC) interactions riNc x N/L 2 , whereas for the disappearance channel we 
neglect backgrounds. Note, that the efficiency e and the cross section a appear as product 
and hence we define an effective cross section 

&u a '■= 0v a *v a ■ (7) 

Many of the quantities — most importantly cross sections and fluxes — appearing in 
Eqs. (JTJ) to ([HI) are subject to (sometimes large) uncertainties. Therefore, the aim is to 
use data from the ND in order to predict the signal in the FD, reducing the dependence on 
external information as much as possible. This can be done efficiently for a disappearance 
measurement. Using Eqs. (pO) and (j3J) one finds 



<° = < D ™^K-^)- (8) 



,FD _ _ND ^FD £ N D 
-And -^FD 

Hence, assuming that the uncertainty on A" fd /A"nd x £ ND /L FD is negligible, a complete 
cancellation of all systematical errors happens (in this idealised discussion), since the same 
combination of a and $ appears in ND and FD. This result is well known and has been 
exploited with success in the K2K [6!] and MINOS [8] experiments. 2 Also, it is the basis 



2 Certainly, in real life the situation is much more complicated than suggested by Eq. (JSJ which should 
illustrate the principle. In the K2K and MINOS analyses many additional complications have been taken 
into account. For example, in both cases the fluxes are rather different in ND and FD, and effects of the 
energy spectrum are included via a near-to-far extrapolation matrix. 
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for all of the latest generation of reactor neutrino experiments like DoubleChooz [19|, [20 . 
For reactor experiments, the idea is that by careful design and construction, the near to far 
comparison can be made so precise that Eq. ([5} holds to an accuracy of better than 1%. 3 

However, the situation is very different for an appearance experiment. Depending on the 
relative importance of the two terms in Eq. (j3J) we can identify two qualitatively different 
regimes for the appearance measurement: first, the regime close to the sensitivity limit of 
the experiment (i.e., small sin 2 26^ ), where the v e events are dominated by background, and 
second, the regime of large sin 2 2#i3, where the actual appearance signal dominates over the 
background in Eq. (j3J). Hence, depending on the regime, one expects that either the error 
on the background or on the signal is most relevant for the sensitivity. As the numerical 
calculations will show, for T2HK the transition between the two regimes occurs roughly at 
sin 2 2^13 ~ 0.01. 

In the background dominated case of small sin 2 2#i3 the ND plays a crucial role in measur- 
ing the background. Under the assumption that the NC background is the same in ND and 
FD, and neglecting the small effect of oscillations on the beam background, P[y e — > i/ e ) « 1, 
the background in the FD can be predicted by the v e events in the ND from Eq. (J2]): 

AT T 2 

< D ' bs = < D ^ff^ (9) 

It becomes important how well the above assumptions are fulfilled, for example, how well 
one can extrapolate the NC background from ND to FD. Moreover, the statistical precision 
for ra^ D is an issue, i.e., the size of the ND, since the beam contains only a small com- 
ponent of z/ e , at the level of 1% of v^. Note that in the discussion leading to Eq. (Q we 
have neglected a background coming from charged- current interactions. This background 
is very different in ND and FD because of oscillations of with sin 2 2 #23 — 1- Therefore, 
this additional background component will further complicate the extrapolation of the back- 
ground measurement from the ND to the FD. While we neglect such a background in the 
current section for simplicity, we do include a background from fi/e miss-identification in 
the numerical calculations, see Sec. HVl 

The signal dominated regime of large sin 2 2# 13 is probably the more interesting case, since 
it will allow for high precision measurements such as CP violation (CPV), and therefore the 
main focus of our work is on this case. From Eq. (j^D one can see that for the signal the 
combination x a Ue is relevant, which cannot be determined by the ND, and Eqs. ([T]) and 
(E} combine to 

FD,sig _ ND ^FD -^ND ^Ve 



«:r — -^ -P- P{v„ -> u e ) . (10) 

Le effective cross sections, namely the ratio 
&v e /&i/u survives. The ability to discover CPV largely depends on the ability to compare the 



iV N D L FD a Vfl 

Obviously, there remains some dependence on the effective cross sections, namely the ratio 



3 For reactor experiments, this very high accuracy is possible because the fiducial mass will be determined 
within less than 0.5%. Moreover, the efficiencies are close to 100%, hence there will be only a very small 
error due to event reconstruction. 
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neutrino and anti-neutrino appearance signals, thus it is useful to look at the ratio of the 
corresponding event rates: 



n FD,sig _ n ND ^ ^ ^ 



FD,sig „ND ^ ^ nZ-n , 77 \ " 

From this discussion we learn that one of the following combinations of quantities has to be 
known in order to predict the signal for the CPV measurement: 

or (<,<) or f^AV (12) 

The first two combinations follow from Eq. ffTT]) : if either the flavour ratio of effective 
cross sections for neutrino and anti-neutrinos separately or the neutrino/anti- neutrino ratio 
for v e and separately are known with good precision then the high statistics and 
77^ samples from the ND allow to predict the CPV signal in the FD. Note that this does 
not require knowledge on the double ratio {b v j 'b v ) / \b v j 'b-p) ■ The last combination in 
Eq. ffl2|) follows directly from Eq. (jSJ): If and fluxes, as well as u e and V e effective 
cross sections are known the signal can directly be predicted without the need of the ND. 
Let us mention again that in Eq. ffl2|) always the product of cross sections times efficiencies 
appears. Uncertainties on both of them contribute to the error on b. Although the preceding 
discussion is highly simplified it captures quite well the behaviour of the near/far detector 
system for an appearance experiment. Many results of our numerical simulation can be 
understood qualitatively with this kind of reasoning, and in the course of our discussion 
we will refer frequently to the arguments presented in this section. In appendix [B] we 
demonstrate that these arguments can be used to approximate the full ND/FD simulation 
by a rather simple (in what concerns the systematics treatment) FD-only setup. 



III. NEUTRINO CROSS SECTIONS 

Since cross section ratios play such an important role in the discussion of systematics, we 



would like to make a few remarks. Based purely on existing experimental data [21[ , without 
the use of a specific model, the errors are in the range 20 — 50%. Especially anti-neutrino 
cross sections are not well measured or for some energy ranges not measured at all. T2HK 
operates in the energy range from 400 — 1200 MeV, therefore most events (at least after the 
single ring cut) will be due to quasi-elastic (QE) reactions and we will focus on these. The 



theory for neutrino scattering off a free nucleon is well understood 22] . Here the cross section 
is given as a function of various form factors 4 , which are all but one well measured. The one 
'free' parameter is the axial mass tua- Based on that formulation, one would expect that the 
ratio of a Ve to <j v can be accurately computed. However, in most detectors the bulk of the 



There are some deviations from pure dipole form, which, however, can be accounted for, see e.g. 23 1 
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fiducial mass stems from heavier nuclei, therefore nuclear effects are non- negligible. Many 
experiments use the Smith-Moniz formalism 24] to account for nuclear effects. Here, the 
nucleus is described as a Fermi gas of nucleons and two new parameters enter: the Fermi 
momentum hp and the binding energy Eb- This description thus introduces a smearing 
of the energy of the outgoing neutrino due to Fermi motion and a reduction of the cross 
section due to Pauli blocking. There are basically two differences between electron and muon 
neutrino scattering: one is the pure kinematic effect of the difference in m e and m M . This 
effect is trivial to account for. The other one is, the fact the momentum transfer to the 
nucleus is going to be different. In order to compute the resulting effect on the ratio it is 
necessary to know the momentum distribution of the bound nucleon. 

Moreover, it is not clear how well the Fermi gas model actually captures the physics of 
the interaction of a neutrino with bound nucleons. Until recently the consensus value for 
the axial mass was tua 
reported by K2K 



26 



1.025 ± 0.021 GeV [25] . This is somewhat in contrast with values 
of m A = 1.20±0.12GeVand MiniBooNE Q of m A = 1.23±0.20 GeV. 
One possible explanation, put forth in 27|], is that the old data was mainly obtained on 
Deuterium, where nuclear effects are very small, whereas the new data used Oxygen (K2K) 
or Carbon (MiniBooNE). This in turn would indicate that there are nuclear effects which are 
not properly included in the Smith-Moniz formalism and thus the value of tu a de termined 
from nuclei is in reality an effective parameter. However, it was pointed out in 23f| that m a 



28 



should be the same or decrease in a nuclear target compared to deuterium, see also 
Thus the experimental situation on the quasi-elastic cross section itself is somewhat unclear. 
It may be that the data on which the value for ttla in 25] is based are not pure QE events, 
which of course would introduce a bias into the determination of of m^. The value of 
itself is not expected to have a large impact on the ratio of cross sections. This example 
is intended to show that there are still open issues in seemingly well understood neutrino 
interaction processes. Also a comparison of several state of art event generator in the range 
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< E v < 2 GeV yields errors in the range 5 %-15 % 

Clearly, from a purely theoretical point of view any correction due to finite lepton masses 
should be small especially for energies E ^ m M , which is the case for T2HK. Surprisingly 
we found only very little literature discussing the cross section ratio. There are many papers 
computing either the z/ M or v e QE cross section on various nuclei, but only in Refs. 



30, 



31 



32 



31 



the error on the ratio is 



we could find a result for both, and v e on Oxygen 5 . In Ref. 
quantified explicitly, however only in the energy range below 500 MeV. Their result is, that 
within the model given in 33] the errors are about 1% coming from the uncertainties of the 



input parameters. They also estimate that physical effects not accounted for introduce no 



more than 5% error on the ratio. Ref. 32j gives both cross sections in the energy range 



relevant for T2HK. The authors of [32] kindly provided their results and we could compute 



5 We do not claim that our survey is complete, but it certainly is representative of the small number of 
pertinent results compared to the total amount of literature about neutrino cross sections. 
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the ratio and compare it to the results in [331] - At 450 MeV we find a difference in the ratios 



33 



of about 3%, which falls within the error estimate given in 

We furthermore compared the results for the ratio obtained with the event generators 
NUANCE Q, GENIE Q and NEUGEN Q and we find between 400-1200 MeV a spread 



of about 1% 6 . With respect to the theoretical calculation in [32] we find a difference of 3% 



ad, 



31 



at 400 MeV which decreases down to 0.5% at 1200 MeV. We extracted the data in 
from the published plots and, within the errors this inevitably introduces, they show a 
similar spread. Summarising, all theoretical sources for the ratio we could find, showed a 
spread of 3% or less throughout the energy range relevant for T2HK. Let us mention that 
in the energy region of a few 100 MeV, the error on the ratio from theory reaches more than 
10%. This energy range is relevant for a beta beam with a Lorentz-7 of around 100 or the 
SPL superbeam 37|. Therefore, in these cases present theory calculations do not provide 
a relevant constraint on the ratio. Also note, that the spread for the ratio of neutrino to 
anti-neutrino cross sections is found to be larger than 10% throughout the whole energy 
range for T2HK. 

MiniBooNE offers an excellent case study of how a recent experiment deals with these 
issues. MiniBooNE is a v e appearance experiment in the same energy region than T2HK, 
which uses, instead of a near detector, the unoscillated sample in the far detector to 
predict their v e backgrounds and also the v e signal. A dedicated investigation of the effect 



of cross section uncertainties has been performed, and the results of Ref. [38| indicate that 
the Vy, to u e cross section ratio has an error of about 8 — 9%. 

The question is, whether one trusts theory calculations, which state that the ratio is 
known to better than 3%. If this is the case, the cross section ratio would have only a small 
impact on the overall error budget and the near-far comparison would effectively control the 
systematics also in an appearance experiment. On the other hand, one should acknowledge 
the fact that neutrino scattering data is sparse and no theory of quasi-elastic scattering 
has been experimentally tested with a large degree of accuracy. Thus there is, at least in 
principle, considerable room for surprises and consensus is needed on whether this risk is 
tolerable in view of a large scale project such as considered here. For instance, there is a 
long standing excess of u e events in sub-GeV atmospheric neutrino data whose origin so 
far is not understood and might reflect a problem in the z/ e /z/ /t cross section ratio. 7 Also, 
the error estimate obtained by MiniBooNE clearly points towards much larger errors of the 
order 10%, than our survey of theory results has found. 

Summarising, it seems that the 10% default errors used here for the individual cross 
sections are somewhat optimistic, especially if one keeps in mind that basically all existing 
data is nearly exclusively for v^. On the other hand, with our defaults and assumption of 
uncorrelated errors for all cross sections, the ratio o Ve jo v ^ would have an effective error of 
~ 10% (see appendix [B]) , which may be larger than expected from theory but is in agreement 



6 NEUGEN and GENIE give nearly identical results. 

7 We thank E. Lisi for pointing this out. 
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with the MiniBooNE numbers. In the following we will take the conservative view point, that 
we would like to be independent of theoretical arguments about the cross section ratio and 
use (by then) available experimental data to control systematical uncertainties. Nonetheless, 
we show results for various constraints on this ratio. Note, that the error relevant for the 
oscillation analysis is the error on the ratio of the product of cross section and efficiency 
Thus, even if there is a tight constraint on the cross section ratios, the efficiencies still will 
need to be determined accurately by other means. 



IV. DESCRIPTION OF THE SIMULATION 

Let us now give some key features of our numerical simulation. More technical details are 
deferred to appendix [A] We consider the following standard setup for the phase I (phase II) 
T2K (T2HK) configuration. The fiducial far detector mass is 22.5 kt (500 kt) and the beam 
power is 0.77 MW (4 MW). The running time is for T2K and T2HK 2 yr for the neutrino 
and 6 yr for the anti-neutrino beams. The baseline is 295 km and we use an average matter 
density of 2.8 gem -3 . This setup is based on [16]. Details of the T2KK setup are given in 
Sec. IV Dl We consider the signals from and v e single ring events, for both the neutrino and 
anti- neutrino run, i.e., disappearance and appearance channels, and we take into account 
also the effect of oscillations on various background components. Energy resolution is treated 
with migration matrices including nuclear effects as well as the contamination of the single 
ring sample with non-quasi-elastic events. We include the intrinsic u e , V e and I7 M (and the CP 
conjugate ones for the anti-neutrino run) as well as neutral current backgrounds. We restrict 
the analysis to reconstructed neutrino energies from 0.4 — 1.2 GeV. This range is divided 
into 8 equidistant bins. For the near detector we assume a water Cerenkov detector with 
fiducial mass of 0.1 kt and otherwise identical properties to the far detector. Specifically, we 
also assume the same acceptance, i.e. a flat near-far ratio. This assumption translates into 
the requirement that the near detector distance must be large enough in order to see the 
decay pipe as point source like the far detector. To be specific, we follow the choice of the 



T2K collaboration and use a baseline of 2 km for the near detector 39( . Further details on 
the detector simulation are to be found in appendix lAl 

The x 2 computation is based on a standard Poissonian form and we use the so-called pull 



approach 40j, |41[ to include the various sources of systematical errors. For the implementa- 



tion of the systematical errors and an explicit definition of the x 2 -function see appendix IA21 



The calculations have been performed with the GLoBES software 42j, |43|, |44j, exploring the 
possibility of user-defined \ 2 m order to include the various systematics with the proper 
correlations. A GLoBES gib-file for our T2HK simulation with an effective systematics 
treatment (see appendix [B]) is available at [44j. 

We include 27 uncorrelated errors listed in Fig.[TJtogether with the adopted default values. 
They include detector normalisations and energy calibration errors, uncertainties on the 
initial fluxes, cross section uncertainties which in our convention include also uncertainties 
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FIG. 1: List of the systematical uncertainties and the adopted default values, as well as the impact of 
various systematics on the T2HK sensitivity to CPV for sin 2 2(?i3 = 0.03. The abscissa shows the smallest 
5qp in [0, it/ 2] for which CPV can be established at 3c We show how the sensitivity is affected if each of the 
27 pulls is switched off (red) or the error is multiplied by 5 (blue). The lower 6 rows show the impact when 
certain combinations of pulls are constrained at 2%: the ratio of v e to cross sections times efficiencies 
(for neutrinos and anti-neutrinos), the ratio of neutrino and anti- neutrino cross sections times efficiencies 
(for e and /x-like events) , the product of flux times v e cross section times v e efficiency (for neutrinos and 
anti- neutrinos), the ratio of e to /i fluxes (for neutrinos and anti-neutrinos), z/ e and F e cross sections times 
efficiencies, and 77 p cross sections times efficiencies. 
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on the efficiencies, as well as errors on NC backgrounds and muon miss-identification. A more 
detailed discussion is given in appendix IA2I including also some motivations for our default 
values. We stress that our defaults should by no means be considered as the most realistic 
values, especially at the time when the experiment is actually performed. Our assumptions 
are motivated by the present situation (see appendix for references) or in other cases are 
very conservative guesses. The purpose of our work is not to advocate specific values for the 
systematics. Instead we want to identify the crucial uncertainties which cannot be eliminated 
with the help of the ND, and hence, for which solid external information is required. 



V. RESULTS 

A. CP violation at T2HK 

In Fig. [2] we show the 3cr sensitivity of T2HK for CPV. We restrict the analysis to the 
range < 5^p c < ir/2, and we neglect the sign(Am| 1 ) degeneracy. When calculating the 
X 2 for 5qp = and 7r we fix all oscillation parameters except from # 13 to their assumed 
true values Am 2 31 = 2.4 x 1(T 3 eV 2 , sin 2 6 23 = 0.5, Am 2 21 = 7.9 x 1(T 5 eV 2 , sin 2 6 12 = 0.3. 
This allows to focus exclusively on the impact of systematics for the discovery of CPV. 8 The 
figure shows the sensitivity from statistical errors only (lower black curve), which is obtained 
by fixing all 27 systematic pulls to zero, as well as the sensitivity for our default choice of 
systematical errors according to Fig. [TJ (upper black curve). Clearly the impact is rather 
dramatic for large 9 13 , in the region sin 2 2# 13 > 10~ 2 . The effect of the various systematics 
in that region is illustrated in Fig. [TJ where we show what happens to the smallest 5qp for 
which CPV can be established at 3a if each single pull is switched off one by one (red bars) 
or multiplied by a factor of 5 (blue bars), assuming sin 2 26 , 13 = 0.03. No single pull has a 
large impact on its own, which highlights the importance of a comprehensive treatment of 
a large number of possible error sources. 

In Fig. [2] we show the case when rather precise information (a hypothetical 1% error) 
is available for either the and z/ M effective cross sections or the v P and V P , effective cross 
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sections. Future cross section experiments such as e.g. MINERVA [45J or SciBooNE [46 
for a 5% accuracy on the absolute cross section. In section IV Bl we explore under which 
circumstances the near detector itself can perform an accurate cross section measurement. 
Note, that the effective cross section is defined as the product of cross section and efficiency. 
Therefore, also the efficiencies would have to be known to better than 1%. Apparently only 
a marginal improvement is possible for precisely known effective cross sections. This is 
not very surprising since for z/ M the near detector can indeed cancel a large fraction of the 
associated errors. Knowing the v e effective cross section would be helpful (c.f. also Fig. [1]), 



For the sensitivity to CPV as shown in Fig. [2] the impact of the uncertainty on the oscillation parameters 
(as well as on the matter density) is negligible. See also the discussion later in connection with Fig. [51 



11 




Q I I I I I I I I I I I I I I I I I I I I I I I I I I 

io" 3 io" 2 10" 1 

sin 2 29 B 

FIG. 2: T2HK CPV sensitivity at 3cr for our default choice of systematical errors according to Fig. [1] 
and for statistical errors only (curves delimiting the shaded region). We show also the sensitivity if certain 
constraints on the product of cross sections times efficiencies a are available: 1% accuracies on a Ve and 
a Vfx for neutrinos and anti-neutrinos, and 5%, 2%, 1% accuracies on the ratios d v ^/a v for neutrinos and 
anti-neutrinos. 

since here the near detector provides very limited information, and the signal is directly 
proportional to this cross section. However, clearly this information alone cannot resolve 
the bulk of the systematics problem. 

In contrast, the situation improves significantly if external information is available on 
the ratio of the effective cross sections gr = d Va ja Vtl for neutrinos and anti- neutrinos. Such 
information can come either from theoretical calculations (for the cross section only, see 
Sec. IHIj) or dedicated experiments. In order to perform this analysis we take into account 
a correlation matrix between the pulls corresponding to the effective cross sections, which 
imposes a constraint on the ratio. Hence, we replace the uncorrelated penalty terms for 
the pulls in Eq. (IA3I) by a matrix correlating the relevant pulls. In Fig. [2] curves are shown 
corresponding to an error on or of 5%, 2% and 1%. Clearly, at the largest values of #13 
the error budget is dominated by <tr, and constraining this quantity basically would allow 
to recover most of the statistical accuracy of the experiment. This is in agreement with the 
discussion in Sec. [Hi and can be understood by the fact the relative size of CP effects is 
smallest for large #13 and hence the absolute accuracy of the prediction of the number of 
oscillated v e events becomes very important. Any CP effect has to be uncovered from this 
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FIG. 3: CPV sensitivity at 3tr as a function of exposure for a true value sin 2 26*i3 = 0.03 for our default 
choice of systematical errors according to Fig. [T] and for statistical errors only (curves delimiting the shaded 
region). The ratio of neutrino to anti-neutrino running is kept constant at 1 : 3. Furthermore we show 
the sensitivity if certain constraints on the product of cross sections times efficiencies a are available: 1% 
accuracies on 5v e and a Vtl for neutrinos and anti-neutrinos, and 5%, 2%, 1% accuracies on the ratios o Ve ja v 
for neutrinos and anti-neutrinos. 

number or, more precisely, from its error. The contribution of backgrounds to the number of 
u e events is small and the contribution to the error is even smaller. The effect of a constraint 
on &R is also shown in the lower part of Fig. [U where we show the impact of constraining 
certain combinations of systematics. There we display also the impact of a constraint on the 
ratio of effective neutrino and anti-neutrino cross sections, which has a very similar effect as 
a constraint on the flavour ratio (even slightly more effective) for restoring the statistics-only 
sensitivity, in agreement with the discussion related to Eq. (|12|) . 

Fig. [3] shows the luminosity scaling of the sensitivity to CPV in the region of large #13 
(sin 2 2#i3 = 0.03) with the same set of curves as in Fig. [2J For this analysis we scale 
simultaneously the beam power and the FD mass between the T2K (0.77 MW, 22.5 kt) and 
T2HK (4 MW, 500 kt) benchmarks defined in Sec. UV[ and we always assume a neutrino (anti- 
neutrino) running time of 2 (6) years, and hence a factor 8 yr is included in the horizontal 
axis in Fig. [31 The ND mass is fixed at 0.1 kt. Obviously, the impact of systematics becomes 
larger as the luminosity increases. Increasing luminosity is a possibility to compensate for 
systematical errors, though a very costly one. For example, if an external constraint on 
the ratio <J Ve j& v at 2% were available T2HK could have 1/10 of the luminosity and still 
achieve the same sensitivity to CPV as in the default configuration. Thus, any optimised 
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FIG. 4: T2HK CPV sensitivity at 3er for a total rate measurement only (left) and without a near detector 
(right) for our default choice of systematical errors according to Fig. [T] and for statistical errors only. The 
dashed curves correspond to an external accuracy of 1% on the ratios d Ve ju v for neutrinos and anti- 
neutrinos. The shaded regions correspond to our standard analysis and are identical to the one in Fig. [2] 

experimental strategy has to find the right balance between spending on measures to mitigate 
systematical effects and on maximising the total luminosity. This may require a dedicated 
effort since some of the necessary measurements may be external to the actual oscillation 
experiment. 

Let us briefly comment on the region of small # 13 close to the sensitivity limit. Fig. [2] 
shows that in this case a constraint on a Ve /a v cannot reduce the effect of systematics. The 
reason is that in this region the precision on the background determines the sensitivity, see 
the discussion in Sec. [TTJ This leads to rather different requirements than in the region of 
large 13 . The relevant question is here to which accuracy the ND data can be used to 
predict the background in the FD, c.f. Eq. (Q. This ability is limited by backgrounds in the 
ND (e.g., NC, muon miss-identification), as well as by statistics in the ND. 

Fig. H] (left) shows the impact of spectral information. The first observation is that the 
pure statistics sensitivity hardly changes if only rate information is used, which shows that 
the main oscillation physics is captured just in the total number of events. 9 However, the 
spectrum is important to disentangle oscillation effects from systematics, especially in the 
regions of very small #13 (to measure the background in the ND) and very large #13 (to avoid 
confusion of the CPV signal with systematics). This result indicates that the systematics 
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Note that spectral information is crucial for resolving the so-called intrinsic degeneracy, see e.g. [17l.l37l.l40|. 
However, since the intrinsic degeneracy does not confuse CP violating and conserving values of i5cp it does 
not affect the sensitivity to CP violation shown in Fig. |4j 
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question in the context of a wide band beam 47|, |48|, |49( might be different, and it would be 



desirable to have a similar analysis also for such a facility. 

In Fig. H] (right) we show the sensitivity to CPV without any near detector. This plot 
highlights the importance of the ND for small #13, where it is needed to constrain the 
background. However, the impact of the ND is somewhat smaller for large #13, sin 2 2^13 > 
0.05, since here the question of backgrounds is less important, whereas the main uncertainty 
comes from the combinations given in Eq. (fl2|) . e.g. the ratio o v<t ja v ^ for which the ND 
provides only rather poor constraints. 

B. Constraints on neutrino fluxes and properties of the near detector 

Let us now discuss the impact of some external knowledge on the initial neutrino fluxes. 
In this case some of our default values might appear slightly too conservative. Information 
on the fluxes requires careful instrumentation of the beam line and data from dedicated 



Hadron production experiments such as MIPP 50j in the case of MINOS, HARP 



for K2K and MiniBooNE, or NA61/SHINE [53j for T2K. For example, in MINOS the goal 



51 



52 



is to constrain at the 5% level using MIPP data. It is beyond the scope of this work 
to do a detailed study of how well neutrino fluxes can ultimately be determined. Here we 
investigate the case of perfectly well known fluxes and how useful that would be for the CPV 
measurement in T2HK. 

The dashed curve in Fig. [5] corresponds to the case of perfectly known fluxes (including all 
sub-dominant flavours in each beam) with all other systematics at the default values. The 
sensitivity improves slightly for sin 2 26 \% > 0.01, but clearly this information is not enough to 
significantly address the systematics problem. The reason can again be understood from the 
discussion in Sec. HO Eq. (flUj) shows that the uncertainty on o v Ja Vfi remains, irrespectively of 
the uncertainty on the fluxes. As mentioned in the paragraph after Eq. (fl2|) flux information 
is only useful in combination with a constraint on the effective v e cross section. This is 
confirmed by the blue curves in Fig. [5l which show that the impact of systematics can 
efficiently be reduced by accurate external information on both, $ V/1 and a Ue (for neutrinos 
and anti-neutrinos). Note that this information is not provided by the ND, but has to come 
from sources outside the considered setup. This seems especially difficult for v e and V e cross 
sections. 

Let us elaborate more on the somewhat surprising result, that knowing the fluxes has 
such a small impact. Indeed, in this case one may expect that the ND should be able to 
provide a measurement of a Ue via the intrinsic v e component in the beam (which is assumed 
to be known perfectly). This argument is true in principle, however, we find that within 
our implementation the NC background and the muon miss-identification in the ND plus 
statistical errors in the ND are enough to spoil this measurement. We illustrate this effect 
in Fig. [6] by showing the CPV sensitivity for sin 2 2#i3 = 0.03 as a function of the ND mass, 
while keeping our T2HK default values for beam power, FD mass, and vjv running times 
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FIG. 5: T2HK CPV sensitivity at 3c. We shows the impact of perfectly known fluxes, as well as constraints 
at 5%, 2%, 1% on a Ue and fluxes, both for neutrinos and anti-neutrinos. The shaded region corresponds 
to our standard analysis and is identical to the one in Fig. [2l 

constant. First, we note that for all systematics at default the sensitivity depends very little 
on the size of the ND, which is consistent with Fig. H] (right). Second, the curve for known 
fluxes (green/dash-dotted) shows a modest gain in sensitivity at the default ND mass of 
0.1 kt from Fig. 0, and some improvement with increasing the mass. However, the situation 
clearly improves if a "perfect" ND without muon miss-identification and NC background is 
assumed. In this case it only depends on the statistical errors in the ND and on the a priori 
accuracy of the fluxes, how well a Ve can be determined by the ND. This can be seen from 
the blue/solid curves in Fig. El corresponding to a "perfect" ND plus a constraint on the 
fluxes, where the accuracy indicated in the figure is implemented as an uncorrelated error 
on each of the flux components. From the plot we find that for a flux uncertainty of 1% the 
pure statistics sensitivity is nearly reached for ND masses of about 1 kt. Let us add that 
such a "perfect" ND would also improve the sensitivity to CPV at small #13, since it would 
provide an accurate determination of the background. 

To summarise this discussion, if precise information on fluxes is available (including the 
v e and V e components) it is possible, in principle, to measure the electron cross sections in 
the ND. This would lead to a strong reduction of the systematics impact, since then both 
$„ and er„ e were known, which is one of the "magic" combinations identified in Eq. (1121) . 
To achieve this situation, in addition to the flux information a careful design of the ND in 
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FIG. 6: CPV sensitivity at 3a for T2HK as a function of the near detector mass for a true value sin 2 2#i3 = 
0.03 for our default choice of systematical errors according to Fig. [1] and for statistical errors only (curves 
delimiting the shaded region). The red/dashed curve corresponds to a ND with perfect e/fi separation and 
with perfectly known NC background, but all other systematics at default. The green/dash-dotted curve 
corresponds to the standard ND but we assume that all fluxes are perfectly known. For the blue/solid curves 
we assume a ND with perfect e/fi separation and without any NC background plus some knowledge on the 
fluxes according to the labels. 



terms of background rejection and its size is necessary. In this respect we mention that for 
T2K a liquid Argon detector with a fiducial mass of 0.1 kt is foreseen at 2 km, which would 
allow to collect a fairly clean sample of v e CC events |39t ] . 



C. Determination of 13 at T2K and T2HK 

Let us now discuss the impact of systematics on the (9 13 measurement. In Fig. [7] we show 
the smallest value of sin 2 26* 13 which can be distinguished from #13 = as a function of the 
luminosity, assuming two representative values for the CP phase which correspond roughly to 
the best and worst sensitivity. The first observation is that for T2K phase I systematics have 
only a small impact, since this measurement is largely dominated by statistics. Numerically 
we find that the sensitivity of T2K decreases from sin 2 26* 13 = 0.0167 to 0.0172 for <5cp = tt/2, 
and from sin 2 2# 13 = 0.0206 to 0.0214 for <5cp = tt. For T2HK systematics have a non- 
negligible impact on the #i 3 discovery reach. The situation is very similar to CPV for small 
# 13 , and the corresponding discussion in Sec. IV Al largely applies also for the #i 3 discovery: 
for this measurement the background dominates in Eq. (J3J), and hence the uncertainty on the 
background is the most relevant systematics. Its impact is controlled by the ability of the 
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FIG. 7: Sensitivity at 3cr to a non-zero #13 as a function of exposure for Sqp — tt/2 and n for our default 
choice of systematical errors according to Fig. Q] and for statistical errors only (curves delimiting the shaded 
region). The ratio of neutrino to anti-neutrino running is kept constant at 1 : 3. Furthermore, we show 
the sensitivity obtained without uncertainty on the intrinsic beam background (by fixing a Ve and the e-like 
fluxes) and without an uncertainty on the NC background in the far detector. 

ND to predict the background in the FD. We show in Fig. [7] also curves assuming a perfectly 
known u e beam background, and no uncertainty at all on the background (i.e., fixing the 
v e beam contamination as well as the NC background). If the total background is fixed 
the sensitivity is close to the pure statistics case. It is interesting to note that for the two 
examples of 5cp shown in the figure the importance of beam and NC backgrounds is different. 
This is an effect of the spectral shapes of the signal relative to the background, since the 
spectrum of the signal depends on the value of <5cp, and also beam and NC backgrounds 
have rather different shapes. 

Fig. [HJ shows the allowed region in the space of sin 2 26^ and #cp obtained by T2K (left) 
and T2HK (right) for some example points of "true" parameter values. As expected, for 
T2K the impact of systematics is small, though not negligible in this case. Furthermore, we 
show that for T2K the uncertainty on the other oscillation parameters has a sizable impact 
on the allowed region. We have checked that this effect comes entirely from the atmospheric 
parameters Am^ and #23- Apparently the disappearance channel does not provide enough 
accuracy on these parameters to avoid an effect on the #13 determination. For the solar 
parameters the accuracy from present data is sufficient to eliminate any effect on the results 
shown in Fig. [HI 

On the other hand, for T2HK the impact of the correlations with the other oscillation 
parameters is negligible, since the high statistics sample of the disappearance channel pins 
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FIG. 8: Allowed regions in the plane of sin 2 2#i3 and Sep for T2K (left) and T2HK (right) for some 
example choices for the input values marked by stars in the figures. We show the allowed regions for all 
combinations of statistical errors only, systematics according to Fig.[TJ all other oscillation parameters fixed, 
and free (where for the solar parameters we impose present errors). For regions labels "osc par free" we 
allow also for a 5% uncertainty on the matter density, which however has a negligible impact on the results. 
The sign(Am§ 1 ) degeneracy is neglected, and #23 UC — ir/4. 

down the atmospheric parameters with high precision. In contrast systematics are much 
more important. For example, for our test point at large 6*13 (sin 2 26^ = 0.03 and 5qp = l.ln) 
the errors on sin 2 2#i3 and Sep are roughly a factor three larger if systematics are included. 
Clearly in this case the inclusion of systematics would spoil the CPV discovery. 



D. T2KK 

We have tested also the case when part of the HK detector is moved to Korea, at a 
baseline of 1050 km. For this analysis we assume that the second FD is located at the 



same off-axis angle as the first one, like in [17|, [54J . For our standard scenario ( "T2KK" ) we 
assume a 250 kt detector, both in Kamioka and Korea. The \ 2 construction for the second 
FD is completely analogous to the ones for the first FD, and we take into account the proper 
correlations of systematics in the three detector system of ND, FD1, FD2. We focus here 
on the CPV discovery in order to compare the T2KK and T2HK performances. Needless to 
say, that a main motivation for having a detector in Korea is to measure the neutrino mass 
hierarchy which we do not discuss here. The hierarchy determination in turn reduces the 
impact of degeneracies, which we have not included here in order to focus on the impact of 
systematics. Such considerations have to be taken into account when evaluating the overall 
potential of T2KK (which is not the purpose of our discussion). 

In the left hand panel of Fig. M the CPV discovery reach of T2KK (blue lines) is shown 



19 



0.5 



0.4 



0.3 



0.2 



0.1 - 




statistics only 



CPV at 3o 



GLoBES 2007 



i i i i 



500 400 
0.25 



0.2 



0.15 



^ CPV at 3a 
\ 

\ 



0.1 - 



0.05 



T 



mass [kt] @ 1050 km 
300 200 100 

r 



T 



T 



rl 
K 



all syst :matics @ default 

X 



10"' 



10" J 




sin 29, 



100 200 300 400 
mass [kt] @ 295 km 



500 



FIG. 9: Left hand panel: sensitivity to CPV at 3tr for T2HK and T2KK. Right hand panel: sensitivity to 
CPV at 3a for two values of sin 2 29i 3 by changing the detector mass between Kamioka and Korea. 



in comparison to T2HK (shaded region). We find that the pure statistics sensitivity is 
slightly worse for T2KK. If all systematics are put at our default values splitting the detector 
yields a somewhat better robustness with respect to systematical uncertainties at large 
sin 2 26* 13 > 10~ 2 . If precise information on cr Ue /a u ^ is available T2HK and T2KK perform 
rather similar. We do not observe a particular cancellation of systematics beyond that 
already present between near and far detector. Having two baselines makes the physics 
signal more distinct and unique and hence it is harder for systematical effects to mimic 
it for large 6* 13 . A similar result was found with a simplified analysis in js5| . Based on 
these results, it seems not necessary to demand that the two detectors are identical or are 
located at the same off-axis angle. This is especially important in the context of the results 
which indicate that a more on-axis location for the detector in Korea would 



m 
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greatly enhance the sensitivity to the mass hierarchy. The larger background present at a 
more on-axis location maybe tolerable, especially if improved algorithms for tt identification 
are used. This issue has been extensively studied in the context of a wide band beam in the 
US 
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In the right hand panel of Fig. we show how the discovery reach changes for various 
distributions of 500 kt fiducial mass between the Korea and Kamioka sites. In the case of 
statistics only (and neglecting the impact of the hierarchy degeneracy) the conclusion would 
be that CPV is best discovered by putting all mass to Kamioka. This conclusion changes in 
presence of systematical errors and now it depends on 8 13 whether T2KK or T2HK performs 
better. 
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VI. SUMMARY AND DISCUSSION 



We have studied the impact of a large number of possible systematical errors on the 
ability of a superbeam experiment to discover CP violation. As a specific example we chose 
T2HK, however our main results should be applicable to all superbeam experiments using a 
narrow band beam. We implemented a realistic description of the far detector and included 
for the first time a near detector in such a study. The emphasis of this work is not to predict 
the actual sensitivity of a given experiment nor the actual size of systematical errors, but 
to show that under semi-realistic assumptions the effects are large and need to be studied 
in more detail. We find that the cancellation of systematics between near and far detectors 
remains incomplete for the appearance channel, due to a lack of information in the near 
detector on the final state. 

In this respect we have identified two qualitatively different regimes depending on the 
size of #13. For small values, close to the sensitivity limit, the main issue is the uncertainty 
on the background. In this case the performance depends on the ability of the near detector 
to predict the background in the far detector. In the regime of large 13 (sin 2 26 13 > 0.01, 
which is probably the more interesting range for this type of experiments) backgrounds are 
a minor issue and the uncertainty on the signal itself dominates. We find that the impact 
of systematics even with a near detector is rather strong in this regime. For instance, for 
T2HK at sin 2 2^i3 = 0.1 the smallest <5cp for which CPV can be established increases from 
0.05-7T for the statistics only case to 0.247T when systematics are included. 

However, we were able to identify crucial combinations of parameters, which, when well 
constrained (at the level of ^ 2%) can restore the sensitivity nearly to its statistics only 
value, namely 

• the ratios of the effective and u e cross sections d v /a v& for neutrinos and anti- 
neutrinos, or 

• the ratios of the effective cross sections between neutrinos and anti-neutrinos, for v e 
and u^, or 

• the initial flux of and the effective v e cross section, both for neutrinos and anti- 
neutrinos. 

With the effective cross section a we mean here the product between physical cross section 
and detection efficiency. The success of a superbeam experiment in the regime sin 2 26*13 ^ 
0.01 will depend to a significant degree on the information available on these combinations. 

Theoretical cross section calculations indicate that the uncertainty on the ratio a v jo Ve 
might actually be at the level of few percent in the T2K energy range of around 700 MeV. 
However, this result has not been tested experimentally. We stress that this would be a 
crucial input in the analysis of a superbeam experiment which is not based on any data. 
Future cross section experiments such as for example MINERz/A may provide a measurement 
of g v at the 5% level. However, from present perspective it seems difficult to obtain a precise 
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measurement for electron neutrino — and especially for electron anti-neutrino cross sections, 
which are essential for predicting the appearance signal. Maybe the only places where these 
cross sections can be measured in the relevant energy range are beta beams or a neutrino 
factory. Note that the absolute normalisation of the cross sections is needed, which always 
is subject to uncertainties on initial fluxes. Precise information on fluxes may be obtained 
from Hadron production experiments, such as MIPP, HARP or NA61/SHINE. 

Apart from CP violation in T2HK, we find that systematics play a minor role for the #13 
discovery sensitivity in T2K (phase I), since this measurement is dominated by statistical 
errors. For the T2KK setup, where half of the HK detector is moved to Korea, we find 
that the second far detector helps somewhat in reducing the effect of systematics on the 
CP violation sensitivity at large 13 . However, this effect does not come from a cancellation 
of systematics, but from a more robust oscillation signal in the very-far detector. Hence it 
seems not necessary to demand that the two far detectors are identical or placed at the same 
off-axis angle. 

In order to focus on the impact of systemtaics we have neglected the hierarchy degeneracy 
in our study. It is known that for T2HK the presence of the degenerate solution reduces 
the sensitivity to CPV in a certain range of the parameter space. We have checked that our 
conclusion on systematics is not changed due to the degeneracy, since also the fake solution 
is affected in a similar way by systematics as the true one. We stress that for a full evaluation 
of the CPV sensitivity and a comparison to other experimental options (such as e.g. T2KK) 
degeneracies have to be included. 

Before concluding we add here a few thoughts on whether and how our results may be 
extrapolated for other facilities, beyond T2HK. Our results indicate that spectral information 
plays an important role in limiting the effect of systematical uncertainties. This suggests that 
the behaviour of a wide band superbeam will be different. Without a detailed simulation 
it is hard to estimate quantitatively whether the impact of systematics is significantly less 
than in the case of the off-axis configuration considered here, and clearly investigations along 
these lines would be an interesting topic for future work. 

For a beta beam in principle similar considerations apply as in the case of the superbeam, 
however there are some important differences. First, the initial flux of electron neutrinos is 
known to good precision. Second, since the signal here is appearance, the relevant cross 
sections are much easier to measure at a MINERz/A type experiment. Hence, it seems easier 
to constrain the beta beam equivalent of the last combination of quantities listed above, 
namely v e fluxes and cross sections. As already mentioned, a close detector at a beta 
beam would probably be an ideal place to measure the electron cross sections needed for a 
superbeam experiment. 

The optimal facility concerning systematics seems to be a neutrino factory. In this case 
intense fluxes of all four flavours Q Ue , $,7 e , <& v ^ are available at the near detector, and they 
are known with very good precision. Hence, all cross sections can be measured accurately at 
the near detector, which allows to predict the appearance signal in the far detector basically 
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free of systematics on fluxes and cross sections. 10 Nonetheless, it would be useful to actually 
prove this by explicit calculation. 
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APPENDIX A: EXPERIMENT SIMULATION AND SYSTEMATICS TREAT- 
MENT 

1. Detector simulation 

The v e and V e appearance signals in the far detector (FD) are calculated in the following 
way. We take into account quasi-elastic (QE) as well as non-quasi-elastic (NQE) charged 



current events using cross sections from the NUANCE v3r503 event generator [34| . However, 
we require that only a single ring is visible in the detector which strongly reduces the number 
of NQE events. At the generator level the requirement is that just one particle momentum 
is above the Cerenkov production threshold. We take into account an energy dependent 
efficiency for e-like events [16]. Since this efficiency is the product of single ring events 
plus particle identification for u e events, we disentangle the one ring efficiency computed 
with NUANCE, properly weighted between QE and NQE events, and extract a particle 
identification efficiency that we assume to be the same for QE, NQE, v e and V e events. The 
absolute efficiency is normalised in order to reproduce the total number of signal events 



provided in Tab. 2 of Ref. [16| for given oscillation parameters and 5 yr neutrino data in 
T2K. We assume the same efficiency function also for anti-neutrino data. For /i-like events 
we take a flat efficiency of 0.9. 

The following background sources are included for the v e appearance signal (and CP 



10 We note that in case of the neutrino factory another important systematics (at large #13) is the uncertainty 
on the matter density Its effect on the CP violation sensitivity has been discussed in Ref. [57( together 
with possibilities to reduce it. 
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symmetric for the V e signal): miss-identified neutral current (NC) events, the intrinsic v e 
and V e beam contamination, V e events from oscillations of the beam component, and a 
tiny background from miss-identified muons from charged current (CC) events (at a rate 
of 0.1%). For the NC background we extract the spectral shape and number of events from 
Ref. [16j by scaling to our exposures. Lacking any detailed information for the anti-neutrino 
mode, we assume the same size and spectrum for the NC background as in the neutrino 
mode. 

For the spectral analysis we take into account the energy reconstruction for QE and NQE 



events via migration matrices calculated by using NUANCE 34J and SK reconstruction 



algorithms (see Ref. [58], p. 139). We use 50 bins in true neutrino energy from zero to 
2 GeV mapped onto 8 bins in reconstructed neutrino energy from 0.4 to 1.2 GeV. In total 
we apply 8 migration matrices: for QE and NQE events for each neutrino flavour z/ e , z/ e , u^, 
V^. Each matrix is normalised to take into account the single ring efficiency The migration 
is consistently applied to signal and v e background events. Precise information on our FD 
simulation including the migration matrices, backgrounds and efficiencies can be recovered 
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from the GLoBES gib-file available at 

For the near detector (ND) we assume the idealised situation that the flux is identical to 
the one of the far detector for no oscillations (i.e., perfect near-to-far extrapolation). For 
definiteness we take a 0.1 kt detector at a distance of 2 km, for T2K as well as for T2HK. 
We use the same migration matrices and efficiencies as in the far detector. For the //-like 
events we assume no background beyond the events from the z/ M and beam components. 
For e-like events we take into account the beam intrinsic v e and V e fluxes, NC events, as well 
as miss-identified muons from (or I7 M ) CC interactions with a rate of 0.1%. For the NC 
events we assume the same spectrum as in the FD with the normalisation scaled according 
to the different mass and baseline of the ND. 



2. x 2 definition and systematics 

For the statistical analysis we adopt the following x 2 function based on Poisson statistics 
in each bin: 



X d . 



ata 



A=l i=l 



7f(0,U-A A + A A ln 



Di 



(Al) 



where the index A runs over the 8 data samples in our analysis obtained by all combinations 
of FD/ND, z//z/-beam, and e//x-like events as given in Tab. [H The samples 1, 2 (3, 4) 
correspond to the appearance (disappearance) channels in the FD, whereas samples 5-8 are 
the ND data. For the T2KK analysis we add 4 more data samples to the ones given in 
Tab. [I] corresponding to the second FD. In Eq. (1A1I) . T^(0,£ a ) is the theoretical prediction 
for energy bin i in data sample A, depending on the oscillation parameters 6 and the pulls 
£ a parameterising the systematic uncertainties. Taking only the leading term in the (small) 



24 



A 


Detector Beam Flavour 


A 


Detector Beam Flavour 


1 


FD 


V 


e-like 


5 


ND 


V 


e-like 


2 


FD 


V 


e-like 


6 


ND 


V 


e-like 


3 


FD 


V 


//-like 


7 


ND 


V 


fi-like 


4 


FD 


V 


fi-like 


8 


ND 


V 


//-like 



TABLE I: Data samples. 



pulls, in general these predictions can be written as 11 

TMQ = Nf{0) • (A2) 

a 

Note that the T t A for the ND (A = 5, . . . , 8) do not depend on the oscillation parameters, 
and hence, the corresponding terms in the x 2 serve only to constrain the pulls £ a . As usual, 
the corresponding "data" Df is taken as the prediction Tf at some assumed "true values" 
for the oscillation parameters, true , and for zero pulls, i.e., Df = Nf-(O tIue ). The final \ 2 
is obtained by adding the penalty terms for the pulls and minimising with respect to them: 

The pull minimisation is performed by using a routine developed by Michele Maltoni. 

In our analysis we include 27 independent pulls to account for systematic uncertainties as 
listed in Fig. [T] together with the adopted default value for the errors a a . They are coupled 
to the theoretical predictions via the couplings 7r^ according to Eq. (1A2|) accounting for 
the correct correlations. The pulls 1 and 2 describe the normalisation uncertainties of FD 
and ND (fiducial mass), correlated between v- and z/-beams but uncorrelated between the 
two detectors. Pulls 3-6 take into account the energy calibration uncertainty, correlated 
between v- and z/-beams but uncorrelated between the two detectors and e- and /i-like 
events. Following Ref. 0] we adopt a value of 2.5%. The pulls 7-16 account for the flux 
uncertainties, which are correlated between the two detectors. For each beam we assume all 
four flavour components to be uncorrelated. For the dominating flux we include in addition 
a linear tilt on the spectral shape. These errors are representative of the situation without 
a dedicated hadron production experiment. They are about the same as K2K would have 
had without the HARP data @. 

Pulls 17-21 parametrise cross section uncertainties. We include an uncorrelated normal- 
isation uncertainty on the total CC cross section for all four neutrino flavours v e , v^, V e , V^. 
Note, that pulls 17-21 also account for the effect of the error on the efficiency. Since we 



In our code we use a slightly more complicated pull dependence in order to make sure that the T t A stay 
always positive, which, however, is equivalent to Eq. (|A2[) at first order in £ a . 



Xdata 



(A3) 
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assume identical detectors for the FD and ND it seems justified to consider the efficiencies 
correlated between the detectors, and hence we consider the pulls 17-20 as the effective 
uncertainty including cross section as well as efficiency uncertainties for the corresponding 
event types. Pull 21 accounts for the uncertainty on the ratio of QE and NQE cross sec- 
tions, which we take fully correlated between flavours and neutrino/anti- neutrino. K2K has 
measured this ratio in different near detectors and has found a spread of 20% among the 
measurements jfij], which we adopt as our default value for the uncertainty. 

The uncertainty of NC events are taken to be completely uncorrelated between ND and 
FD. For the FD we include pull 22 for the NC normalisation correlated between v- and 
77-beams, whereas pull 23 accounts for the relative uncertainty. Pulls 24 and 25 account for 
the NC backgrounds to e-like events in the ND, uncorrelated between v- and z/-beams. And 
finally, we include an uncertainty for the rate of muon miss-identification in the ND, again 
uncorrelated between v- and z/-beams (pulls 26 and 27). For the T2KK analysis we add 3 
more pulls, one for the normalisation and two for the energy calibration of the detector in 
Korea. 

Let us give one explicit example, how the theoretical predictions according to Eq. (1A2j) 
are constructed, e.g., for the v e appearance channel in the FD (A = 1). In this case we have 

= n l^L e + n u^-w e + n NC + n u e -bea,m + n V e -bea,m + n V^V e + n miss-ID • (A4) 

The n l x correspond to the oscillation signal and the various backgrounds as described above. 
For simplicity we suppress the dependence on the oscillation parameters, however, in the 
calculation also oscillations of backgrounds are properly included. Then Eq. flA2[) reads 

T i = Nii 1 + 6) + 6) ^calibFD.e 

+ £7 + n^l + n l NC + n miss -m) + £ 8 Trjtftt 

+ ^9 n l e -bea.m + 60 ^-beam + 6l n v^V e 

+ 67 + + <-beam) + 6s (4 e -beam + < M -P e ) + &9 <iss-ID 

+ 6i(<^-<^ E J/2 

+ 6 2 <c + 63< c /2- (A5) 

Here 7if calibFD e and Ttl tQt account for the energy calibration and z/^-flux tilt, respectively. The 
6 lines take into account normalisation and calibration, z/^-flux uncertainty, uncertainties of 
the other flux components, total cross section and efficiency uncertainties, the QE/NQE 
ratio, and NC uncertainties. The last line for the T7-beam (A = 2) would read +£22 n l NC — 
63 n Nc/2 hr order to account for the correlations of the NC pulls as described above. The TV 4 
for the other samples are defined in an analogous way. Through this construction we make 
sure that we use only the information which is actually provided by the ND measurements. 
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io" 3 io" 2 10" 1 

sin 2 29 13 

FIG. 10: Comparison of the effective FD description (dashed curves) with the full ND/FD setup. The 
shaded region, as well as the blue curves correspond to our standard simulation and are identical to the 
corresponding curves in Fig. O 



APPENDIX B: A FAR DETECTOR-ONLY SETUP WITH EFFECTIVE SYS- 
TEMATICS 



Previous sensitivity studies (such as for example the ones in Refs. (la, |37|, |40J) do not 
include a ND in the simulation, and some values for systematical uncertainties are adopted, 
which are assumed to implicitly encode information from the ND. In this appendix we ex- 
amine in which way such choices for systematical errors should be interpreted. We calculate 
the sensitivity to CPV for T2HK for a FD-only configuration (denoted by FD efr ), by using 
exactly the same detector simulation and backgrounds as before, but we include only four 
independent effective systematical uncertainties: the normalisations of the appearance sig- 
nal and the normalisations of the total background, both for neutrinos and anti-neutrinos 
(<Ty lg , cr^ lg , cx^ g , <% g ). According to the discussion in Sec.QTJone expects that a^ s and a^ g will 
be relevant for the sensitivity at small #13, whereas a^ e and a^ g will dominate the sensitivity 
at large #13, and this is indeed the behaviour we find. 

In Fig. [10] we compare this FD eff simulation with the full ND/FD configuration. As 
it must be, the pure statistics sensitivities are identical. It turns out that the ND/FD 
sensitivity with all systematics at the default values according to Fig. [TJ is reproduced rather 
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accurately by FD for the following choice of systematics: 

a si g = ^ig = 1Q% ; a b g = a H = 3 5% (B1) 

We conclude from these numbers that our specific implementation of the ND provides a 
3.5% measurement of the background, whereas the effective error on the signal turns out to 
be 10%. To reproduce the curves corresponding to a constraint on the o Vf ,jo Vv ratio at x% 
in the full ND/FD case, one simply has to set af g = cr" s = x% for FD cff . This confirms 
the arguments given in Sec. [Til that the error on the ratio a Ue /a Uv directly translates into an 
error on the appearance signal. 

The FD eff calculations can be performed with a standard GLoBES analysis. A gib-file for 
this effective FD simulation for T2HK is available at [441 ] . 
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